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We develop the kinematics in Matrix Gravity, which is a modified the- 
ory of gravity obtained by a non-commutative deformation of General 
Relativity. In this model the usual interpretation of gravity as Rieman- 
nian geometry is replaced by a new kind of geometry, which is equivalent 
to a collection of Finsler geometries with several Finsler metrics depend- 
ing both on the position and on the velocity. As a result the Riemannian 
geodesic flow is replaced by a collection of Finsler flows. This naturally 
leads to a model in which a particle is described by several mass parame- 
ters. If these mass parameters are different then the equivalence principle 
is violated. In the non-relativistic limit this also leads to corrections to the 
Newton's gravitational potential. We find the first and second order correc- 
tions to the usual Riemannian geodesic flow and evaluate the anomalous 
nongeodesic acceleration in a particular case of static spherically symmet- 
ric background. 



1 



1 Introduction 

Gravity is one of the most universal physical phenomenon. It is this universality 
that leads to a successful geometric interpretation of gravity in terms of Rieman- 
nian geometry in General Relativity. General Relativity is widely accepted as a 
pretty good approximation to the physical reality at large range of scales. 

We would like to make two points here. First of all, the experimental evidence 
points to the fact that all matter exhibits quantum behavior at microscopic scales. 
Thus, it is generally believed that the classical general relativistic description of 
gravity is inadequate at short distances due to quantum fluctuations. However, de- 
spite the enormous eff"orts to unify gravity and quantum mechanics during the last 
several decades we still do not have a consistent theory of quantum gravity. There 
are, of course, some promising approaches, like the string theory, loop gravity and 
non-commutative geometry. But, at the time, none of them provides a complete 
consistent theory that can be verified by existing or realistic future experiments. 

Secondly, in the last decade or so it became more and more evident that there 
might be a few problems in the classical domain as well. In addition to the old 
problem of gravitational singularities in General Relativity these gravitational 
anomalies include such eff"ects as dark matter, dark energy. Pioneer anomaly, flyby 
anomaly, and others [[8]|. They might signal to new physics not only at the Planck- 
ian scales but at very large (galactic) scales as well. 

This suggests that General Relativity, that works perfectly well at macroscopic 
scales, should be modified (or deformed) both at microscopic and at galactic (or 
cosmological) scales (or, in the language of high energy physics, both in the ultra- 
violet and the infrared). It is very intriguing to imagine that these eff'ects (that is, 
the quantum origin of gravity and gravitational anomalies at large scales) could be 
somehow related. Of course, this modification should be done in such a way that 
at the usual distances the usual General Relativity is recovered. This condition 
puts some constraints (experimental bounds) on the deformation parameters; in 
the case of non-commutative field theory such bounds on the non-commutativity 
parameter were obtained in Q. 

In this paper we investigate the motion of test particles in an extended theory 
of gravity, called Matrix Gravity, proposed in a series of recent papers [Tl|2l|3]|. 
The motivation for such a deformation of General Relativity is explained in detail 
in |l2l. The very basic physical concepts are the notions of event and the space- 
time. An event is a collection of variables that specifies the location of a point in 
space at a certain time. To assign a time to each point in space one needs to place 
clocks at every point (say on a lattice in space) and to synchronize these clocks. 
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Once the position of the clocks is fixed the only way to synchronize the clocks 
is by transmitting the information from a fixed point (say, the origin of the coor- 
dinate system in space) to all other points. This can be done by sending a signal 
through space from one point to another. Therefore, the synchronization proce- 
dure depends on the propagation of the signal through space, and, as a result, on 
the properties of the space it propagates through, in particular, on the presence of 
any physical background fields in space. The propagation of signals is described 
by a wave equation (a hyperbolic partial differential equation of second order). 
Therefore, the propagation of a signal depends on the matrix of the coeflicients (a 
symmetric 2-tensor) g'^^(x) of the second derivatives in the wave equation which 
must be non-degenerate and have the signature ( — h . . . -I-). This matrix can be 
interpreted as a pseudo-Riemannian metric, which defines the geodesic flow, the 
curvature and, finally, the Einstein equations of General Relativity (for more de- 
tails, see [2J). 

The picture described above applies to the propagation of light, which is de- 
scribed by a single wave equation. However, now we know that at microscopic 
scales there are other fields that could be used to transmit a signal. In particular, 
the propagation of a multiplet of A'^ gauge fields is described not by a single wave 
equation but by a hyperbolic system of second order partial differential equations. 
The coefficients at the second derivatives of such a system are not given by just 
a 2-tensor like g''^ix) but hy a. N x N matrix-valued symmetric 2-tensor a^'^{x). 
If a^" does not factorize as a^'' Eg^^", where E is some non-degenerate matrix, 
then there is no geometric interpretation of this hyperbolic system in terms of a 
single Riemannian metric. Instead, we obtain a new kind of geometry that we call 
Matrix Geometry, which is equivalent to a collection of Finsler geometries. In 
this theory, instead of a single Riemannian geodesic flow, there is a system of N 
Finsler geodesic flows. Moreover, a gravitating particle is described not by one 
mass parameter but by A^^ mass parameters (which could be different). Note that 
because the tensor a^^ is matrix-valued, various components of this tensor do not 
commute, that is, {a^", a"!^] 4^ 0. In this sense such geometry may be also called 
non-commutative Riemannian geometry. In the commutative limit, a^^ g'^'^ and 
we recover the standard Riemannian geometry with all its ingredients. Only the 
total mass of a gravitating particle is observed. For more details and discussions 
see [2ii3J. 

As we outlined above. Matrix Gravity is a non-commutative modification of 
the standard General Relativity in which the metric tensor g^'^ is replaced by a 
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Hemiitian N x N matrix- valued symmetric two-tensor 

a''" = g^^'I + Kh''^ (1.1) 

where I is the identity matrix, h^^ is a Hermitian matrix- valued traceless symmetric 
tensor, i.e. 

gi^y = l.tra"', trh^' = 0, (1.2) 

and ;>i is a deformation parameter. 

The dynamics of the tensor field a^^ is described by a diff'eomorphism invariant 
action, 

S(a) = ^ dx £.{a,da) , (1.3) 

where dx is the standard Lebesgue measure on the spacetime manifold and X(a, da) 
is the Lagragian density. Of course, as % ^ this action should reproduce the 
usual Einstein-Hilbert action functional 



5(g) = J dxg'^\R-2K), (1.4) 

where is a pseudo-Riemannian metric, g = \ detg^y|, R is the scalar curvature 
of the metric g, G is the Newton constant and A is the cosmological constant. 

The action of matrix gravity can be constructed in two different ways. One 
approach, developed in [P, "21, is to try to extend all standard concepts of differ- 
ential geometry to the non-commutative setting and to construct a matrix-valued 
connection and a matrix- valued curvature. 

The second approach, developed in IQ, is based on the spectral asymptotics 
of a self-adjoint elliptic partial differential operator L of second order with a pos- 
itive definite leading symbol <Ti{x,^) = a'^^(x)^^^y. It is well known that there is 
an asymptotic expansion as ? ^ of the L^-trace of the heat semigroup of the 
operator L 

oo 

Tr L2 expi-tL) ~ (4mr'^ t'Aj, , (1.5) 

k=0 

where Ai, are spectral invariants of the operator L. For the Laplace-Beltrami op- 
erator L = -g^^'^fNv these coefficients are well known, and, it turns out that the 
Einstein-Hilbert action is nothing but a linear combination of the first two coeffi- 
cients, that is, 

5(g) = Y^(6Ai-2AAo). (1.6) 



a/2 
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Therefore, a similar functional, which is automatically diffeomorphism-invariant, 
can be constructed by computing the same heat kernel coefficients for a more 
general partial differential operator L of non-Laplace type (for more details, see 

m). 

The field equations for the tensor d'^, that we call non-commutative Einstein 
equations are obtained by varying the action with respect to a^''; in the vacuum we 
have, 

d£ d£ 

where a"^,^ = d^^a"^. 

The action has an additional new global gauge symmetry 

a^'^x)^ Ua^\x)U-\ (1.8) 

where t/ is a constant unitary matrix (for more details, see the papers cited above). 
By the Noether theorem this symmetry leads to the conserved currents (vector 
densities) 



In other words, this suggests the existence of new physical charges 



Q = j dxj^, 

where dx denotes the integration over the space coordinates only. These charges 
have purely noncommutative origin and vanish in the commutative limit. 

One can easily localize this global symmetry by introducing the local gauge 
transformations 

ai'^x) ^ U{x)a''''(x)U-\x) , (1.9) 

where U{x) = expaj(.x;), and oj is an anti-Hermitian matrix- valued function of 
compact support, i.e. vanishing at infinity, and a new Yang-Mills field that 
transforms as one-form under diff"eomorphisms and as a connection under the 
gauge transformation. All the geometric structures, including the connection coef- 
ficients, the curvature etc, become covariant under the local gauge transformations 
if one simply replaces the partial derivatives by d^+S^. This leads to a gauged 
version of the above functionals (for more details see [HI 121 121). 

This model may be viewed as a "noncommutative deformation" of Einstein 
gravity (coupled to a Yang-Mills model in the gauged version), which describes, 
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in the weak deformation limit, a.sK ^ 0, General Relativity, Yang-Mills fields (in 
the gauged version), and a multiplet of self-interacting massive two-tensor fields 
of spin 2 that interact also with gravity and the Yang-Mills fields. One should 
stress here that x is a formal parameter that does not have a particular physical 
value; it is just a tool to develop the perturbation theory in h^'^. At the end of the 
derivation we can just set % = 1 . What is measured and describes the extent of the 
non-commutative deformation is the tensor h'^^(x) and its derivatives, which can 
be parametrized by he invariants of this tensor like g^agv/sjfi^h^'^h"^- 

Our approach should be contrasted with the non-commutative extensions of 
gravity on non-commutative spaces with non-commutative coordinates 

[x^,x''] = e^\ (1.10) 

where is a constant anti-symmetric matrix, and the Moyal product 



fix) ★ g(x) = exp I ^^"^ ^ I /(^ + y)8(^ + z) 



y=z=o 



(1.11) 



This approach immediately rasies the question on the nature of the coordinates x^. 
The non-commutativity condition can only be true for some priviledged coordi- 
nates, like Cartesian (or inertial) coordinates. The condition that the matrix 6^^ is 
constant, that is, 

d„e>'' = 0, (1.12) 

is not covariant. It breaks the diffeomorphism invariance of the theory, and, as 
the result, Lorentz invariance. Thus any such theory cannot be diffeomorphism 
invariant. One could try to replace it by the covariant condition 

V^ff"' = 0, (1.13) 

where Vq. are covariant derivatives with respect to some background metric. Then 
the matrix ff^" could be viewed as a covariantly constant antisymmetric 2-tensor. 
However, the integrability conditions for this equation lead to very strong alge- 
braic constraints on the curvature of the metric. Thus, the breaking of the dif- 
feomorphism invariance (and as a result of Lorentz invariance) is an unavoidable 
feature of this approach to non-commutative gravity. Therefore, such theories can 
be ruled out by very restrictive experimental bounds on the non-commutativity 
parameter [|51. 

By contrast, the status of diffeomorphism invariance (and the Lorentz invari- 
ance) in Matrix Gravity is exactly the same as in General Relativity, namely, both 
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theories are difFeomorphism invariant, so, there are no preferred coordinates, and 
a condition like (|1.10l) is impossible. In our model it is not the coordinates that 
do not commute, but the metric! Therefore, the recent strong experimental con- 
straints on the violation of Lorentz invariance do not apply to Matrix Gravity. It is 
rather the violation of the Equivalence Principle that is critical for Matrix Gravity. 
This feature could be used for an experimental test of the theory in the future. 

One should also mention the relation of our approach to so called "analog 
models of gravity". In particular, the analysis in is surprisingly similar to the 
analysis of our papers [1, 2J. The authors of Pl consider a hyperbolic system 
of second order partial differential equations, the corresponding Hamilton-Jacobi 
equations and the Hamiltonian system as we did in |[Il|2l. In fact, their f^" is equiv- 
alent to our matrix-valued tensor a'^^. However, their goal was very different — they 
impose the commutativity conditions on (eq. (44)) to enforce a unique effective 
metric for the compatibility with the Equivalence Principle. They barely mention 
the general geometric interpretation in terms of Finsler geometries as it "does not 
seem to be immediately relevant for either particle physics or gravitation" The 
motivation of the authors of (H is also very different from our approach. Their 
idea is that gravity is not fundamental so that the effective metric simply reflects 
the properties of an underlying physics (such as fluid meachnics and condensed 
matter theory). They just need to have enough fields to be able to parametrize an 
arbitrary effective metric. In our approach, the matrix- valued field a'^'' is funda- 
mental; it is: i) non-commutative and ii) dynamical. 

The main goal of the present paper is to investigate the motion of test particles 
in a simple model of matrix gravity and study the non-geodesic corrections to 
general relativity. 

The outline of this work is as follows. In Sect. 2. we develop the kinemat- 
ics in Matrix Gravity. In Sect. 3. we compute the first and second order non- 
commutative corrections to the usual Riemannian geodesic flow. In Sect. 4 we 
find a static spherically symmetric solution of the dynamical equations of Matrix 
Gravity in a particular case of commutative 2x2 matrices. In Sect. 5 we evaluate 
the anomalous acceleration of test particles in this background. In Sect. 6 we 
discuss our results. 
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2 Kinematics in Matrix Gravity 
2.1 Riemannian Geometry 

Let us recall how the geodesic motion appears in General Relativity, that is, in 
Riemannian geometry (for more details, see O). First of all, let 

F{x,0=^p^\ (2.1) 

where is a non-vanishing cotangent vector at the point x, and = g>^''(x)^^^y 
(recall that the signature of our metric is (- + . . . +)). Obviously, this is a homo- 
geneous function of f of degree 1, that is, 

F(x,^) = AF(x,0. (2.2) 

Let 

H(x,0 = -\FHx,^) = ^\e'. (2.3) 

This is, of course, a homogeneous polynomial of of order 2, and, therefore, the 
Riemannian metric can be recovered by 

g''\x) = -f—Hix,^). (2.4) 

Now, let us consider a Hamiltonian system with the Hamiltonian H{x, ^) 
dx" dH(x,^) 



dt d^f 



/Wv, (2.5) 



The trajectories of this Hamiltonian system are, then, nothing but the geodesies of 
the metric g^y. Of course, the Hamiltonian is conserved, that is, 

g''\xm,{t)m = -E , (2.7) 

where £ is a constant parameter. 
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2.2 Finsler Geometry 

As it is explained in [Si Matrix Gravity is closely related to Finsler geometry 
||9l rather than Riemannian geometry. In this section we follow the description 
of Finsler geometry outlined in [91. To avoid confusion we should note that we 
present it in a slightly modified equivalent form, namely, we start with the Finsler 
function in the cotangent bundle rather than in the tangent bundle. 

Finsler geometry is defined by a Finsler function F{x, ^) which is a homoge- 
neous function of of degree 1 and the Hamiltonian 

H{x,0 = -\f\x,^). (2.8) 
Such Hamiltonian is still a homogeneous function of of degree 2, that is, 

^,—H(x,0 = 2Hix,0, (2.9) 

but not necessarily a polynomial in ! 
Now, we define a tangent vector u by 

u>' = —Hix,0, (2.10) 

and the Finsler metric ^ 

G'\x,^) = -^Hix,^). (2.11) 

The difference with the Riemannian metric is, obviously, that the Finsler met- 
ric does depend on more precisely, it is a homogeneous function of of degree 
0, i.e. 

G'\x,A^) = G^''(x,^), (2.12) 

so that it depends only on the direction of the covector ^ but not on its magnitude. 
This leads to a number of useful identities, in particular. 



Hix,^) = ^Cr\x,^)^,^y, (2.13) 



and 



= G'\X,^)^y. 



(2.14) 
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Now, we can solve this equation for treating as independent variables to 

get 

^f, = G^y(x,u)u\ (2.15) 
where G^y is the inverse Finsler metric defined by 

G^Ax,u)G''"(x,^) = S"^. (2.16) 

By using the results obtained above we can express the Hamiltonian H in terms 
of the vector m^, more precisely we have 



H(x, ^{x, u)) = ^Gf,yix, u)u^u'' . (2. 17) 



The derivatives of the Finsler metric obviously satisfy the identities 

^G^'(x,0 = ^G^^ix,^) = ^G"^x,^) , (2.18) 

O^a O^p d^y 

^,^G-(x,^) = ^^—G^%x,^) = 0, (2.19) 

and, more generally. 

This means, in particular, that the following relations hold 



O^cc out' 



— = &'"(x, 0, ^ = G^aix, u) . (2.21) 



It is easy to see that the metric G^yix, u) is a homogeneous function of u of 
degree 0, that is, 

u''—GUx,u) = 0, (2.22) 

and, therefore, H(x, ^(x, u)) is a homogeneous function of u of degree 2. This 
leads to the identities 

1 d 

^, = -—H(x,^(x,u)), (2.23) 

2 out' 

1 

Gf,y(x, u) = - ^ H{x, ^{x, u)) . (2.24) 
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Finally, this enables one to define the Finsler interval 



'liv 

so that 



ds^ = Guvix, x)dxf^dx^ , (2.25) 



dr = y-ds^ = ^-G^y{x, x)xi'x' dt = F(x, ^(x, x))dt , (2.26) 

where 



= ^, = G,y(x,x)x\ (2.27) 



dxf" 
dt 

By treating H{x, ^) as a Hamiltonian we obtain a system of first order ordinary 
differential equations 



d^ dH{x,0 
dt d^y, 

d^, dHix,^) 



dt dx^' 



(2.29) 



The trajectories of this Hamiltonian system naturally replace the geodesies in Rie- 
mannian geometry. Again, as in the Riemannian case, the Hamiltonian is con- 
served along the integral trajectories 

H(x(t),m) = -E. (2.30) 

Of course, in the particular case, when the Hamiltonian is equal to H{x, ^) = 

all the constructions derived above reduce to the standard structure of Riemannian 

geometry. 



2.3 Matrix Gravity 

The kinematics in Matrix Gravity is defined as follows. In complete analogy with 
the above discussion we consider the matrix 



A{x,0 = a''\x)^,^,, (2.31) 

where a''" is the matrix- valued metric (11.11) . As we mentioned in the introdution 
this expression has been already encountered in physics, in particular, in [4] it is 
shown that it is the most general structure describing "analog models" for gravity. 

This is a Hermitian matrix, so it has real eigenvalues /i,(x, ^), z = 1, 2, . . . , A'^. 
We consider a generic case when the eigenvalues are simple. We note that the 
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eigenvalues /?,(jc, ^) are homogeneous functions (but not polynomials!) of ^ of de- 
gree 2. Thus, each one of them, more precisely ^-hi{x, ^), can serve as a Finsler 
function. In other words, we obtain A'^ different Finsler functions, and, therefore, 
different Finsler metrics. Thus, quite naturally, instead of a single Rieman- 
nian metric and a unique Riemannian geodesic flow there appears A'' Finsler met- 
rics and A'^ corresponding flows. In some sense, the noncommutativity leads to a 
"splitting" of a single geodesic to a system of close trajectories. 

Now, to define a unique Finsler metric we need to define a unique Hamiltonian, 
which is a homogeneous function of the momenta of degree 2. It is defined in 
terms of the Finsler function as in (12.81) which is a homogeneous function of the 
momenta of degree 1. To define a unique Finsler function we can proceed as 
follows. Let fii, i = I,. . . ,N, be some dimensionless real parameters such that 



;=i 



(2.32) 



so that there sccq (N - 1) independent parameters. Then we can define the Finsler 
function by 



N 



Fix, ^) = J] l^i ^-hi(x, . 



(2.33) 



i=\ 



Notice that, in the commutative limit, as x ^ and a^^ = g^"!, all eigenvalues 
of the matrix A(x,^) degenerate to the same value, hi{x,^) = and, hence, the 
Finsler function becomes F(x, ^) = VH?F- In this case the Finsler flow degener- 
ates to the usual Riemannian geodesic flow. 

Next, we define the Hamiltonian according to eq. (|2.8I) 



H(x,^) 



1 



/ N 



^^liyj-hi{x,^) 



V i=\ 



^ N 

-J^n^hiix,^)- l^if^jyjhi(x,0hj(x,0. (2.34) 

i=l l<i<j<N 



In a particular case, when all parameters //, are equal, i.e. //, = 1 /N, the Finsler 
function reduces to 



(=1 



(2.35) 
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By using the decomposition of the matrix-valued metric d^^ as in one can 
see that 

UrA{x,0 = \^\\ (2.36) 

and, therefore, 

N 



^Yj^t{x,^) = \^\\ (2.37) 

i=l 

Thus, we conclude that in this particular case 



1 ' 



N 

l<i<j<N 



(2.38) 



It is difficult to give a general physical picture of these models since the Hamil- 
tonian is non-polynomial in the momenta. Hamiltonian systems with homoge- 
neous Hamiltonians have not been studied as thoroughly as the usual systems 
with quadratic Hamiltonians and a potential. 

2.4 Kinematics 

The problem is, now, how to use these mathematical tools to describe the motion 
of physical massive test particles in Matrix Gravity. The motion of a massive 
particle in the gravitational field is determined in General Relativity by the action 
which is proportional to the interval, that, is, 

-S* particle ^ ~ J yj-g^y(x)dx''dx'' = - J m yl-\x\^dt , (2.39) 
P\ h 

where m is the mass of the particle. Pi and are the initial and the final position 
of the particle in the spacetime, f is a parameter, t\ and ?2 are the initial and the final 
values, ^ = ^ and |ip = ^^y(x)i'^i'' . This action is, of course, reparametrization- 
invariant. So, as always, there is a freedom of choosing the parameter t. We can 
always choose the parameter to be the affine parameter such that \x\^ is constant, 
for example, if the parameter is the proper time f = r, then |ip = -1. The Euler- 
Lagrange equations for this functional are, of course, 

^ = -^+T\p{x)x'^xP = Q, (2.40) 
dt dt^ 
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where Y^a/i are the standard ChristofFel symbols of the metric g^y. Of course, the 
equivalence principle holds since these equations do not depend on the mass. 

In Matrix Gravity a particle is described instead of one mass parameter m by 
N different mass parameters 

nii = mfii , (2.41) 



J^m^. (2.42) 



where 

N 

m 

i=l 

The parameters m, describe the "tendency" for a particle to move along the trajec- 
tory determined by the corresponding Hamiltonian hi(x,^). In the commutative 
limit we only observe the total mass m. 

We define the Finsler function F{x,^) and the Hamiltonian H(x,^) as in eqs. 
(12.331) and (12.341) . Then the action for a particle in the gravitational field has the 
form 



' particle 



mF(x,^(x,x)) dt. (2.43) 



Thus, the Finsler function F{x,^{x,x)) (with the covector expressed in terms 
of the tangent vector x^) plays the role of the Lagrangian. To study the role of 
non-commutative corrections, it is convenient to rewrite this action in the form 
that resembles the action in General Relativity. 

5* particle = - J m^six, x) y]-\x\^dt , (2.44) 

h 

with some "effective mass " mef{{x, x) that depends on the location and on the 
velocity of the particle 



hi(x,^(x)) 

m^six, x) = ) ^ Mi — . (2.45) 



This action is again reparametrization-invariant. Therefore, we can choose 
the natural arc-length parameter so that F{x, ^{x, i)) = 1 . Then the equations of 
motion determined by the Euler-Lagrange equations have the same form 

^|^+yVx,i)i'^/ = 0, (2.46) 
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where y^ap{x, x) are the Finsler ChristofFel coefficients defined by the equations 
that look identical to the usual equations but with the Finsler metric instead of the 
Riemannian metric, that is, 

yafi(x,X) = -&'\x,^(x,x))i—Gyfi(x,x) + —Gya(x,x) - —Gap(x,X)\ . 

(2.47) 

To study the role of non-commutative corrections it is convenient to rewrite 
these equations in a covariant form in the Riemannian language. In the commuta- 
tive limit, as X ^ 0, we can expand all our constructions in power series in x so 
that the non-perturbed quantities are the Riemannian ones. In particular, we have 

r^pix, x) = r^^(x) + e^„p(x, x) , (2.48) 

where 9^a/3 are some tensors of order x. Then the equations of motion can be 
written in the form 

— =Al,^^(x,x), (2.49) 

where 

_ = _,rv.),.^ e.5o, 

and 

Al,^^(x,x) = -e\p{x,x)x"xP , (2.51) 
is the anomalous nongeodesic acceleration. 



3 Perturbation Theory 

We see that the motion of test particles in matrix Gravity is quite different from 
that of General Relativity. The most important difference is that particles exhibit 
a non-geodesic motion. In other words, there is no Riemannian metric such that 
particles move along the geodesies of that metric. It is this anomalous acceleration 
that we are going to study in this paper. 

In the commutative limit the action of a particle in Matrix Gravity reduces to 
the action of a particle in General Relativity with the mass m determined by the 
sum of all masses m,. In this paper we consider two different cases. In the first 
case, that we call the nonuniform model, we assume that all mass parameters are 
different, and in the second case, that we call the uniform model, we discuss what 
happens if they are equal to each other. 
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3.1 Nonuniform Model: First Order in k 

So, in this section we study the generic case when the parameters are different. 
As we already mentioned above, in this case the Finsler function F{x, £) is given 
by (I2.33D . By using the decomposition (11.11 ) of the matrix-valued metric a^^ we 
have 

A(x, ^) = d'\x)^^^, = + Kh^'\x)^,^, . (3.52) 
Therefore, the eigenvalues of the matrix A{x, ^) are 

hi{x,0 = \^? + ^^i{x,0, (3.53) 

where Ai{x,^) are the eigenvalues of the matrix h^'^(x)^fj^v. In the first order in x 
we get the Finsler function 

F{x,^) = Vh^|i + 0(x^), (3.54) 



and the Hamiltonian 



H(x,^) = ^\^\^ + x^P(x,0 + 0(x^), 



(3.55) 



where 

N 

P(x,^) = Y,fiMx,0- (3.56) 

!=1 

By using the fact that P(x, ^) is a homogeneous function of ^ of order 2, we find 
the Finsler metric 

O'^x, = g'^x) + M^'^x, ^) + 0(x^) , (3.57) 

and its inverse 

G^,y(jc, u) = ^pv(^) - Miivi^-' ^(-^' ")) + Oix}) , (3.58) 

where ^ 

q^\x,0 = l--^P(x,^). (3.59) 

Here the indices are raised and lowered with the Riemannian metric, and 

u^ix, = G^"(x, ^)^y , ^^(x, u) = G^Ax, u)u' . (3.60) 
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Since P(x, ^) is a homogeneous function of ^ of order 2 we have 

P(x,0 = q>'\x,^)^,^,. (3.61) 

Note that since tr/z'^^ = the matrix W^^^^^y is traceless, which implies that the 
sum of its eigenvalues is equal to zero. Thus, in the uniform case, when all mass 
parameters /i, are the same, the function P{x,^) vanishes. In this case the effects 
of non-commutativity are of the second order in x; we study this case in the next 
section. 

We also note that 

I^P = Iwp - 2xPix, ^{x, u)) + 0{x^) . (3.62) 
Thus, our Lagrangian is 

F(x,^(x,x)) = V^(l-%^^^^^^|^j + 0(x2), (3.63) 

Finally, we compute the Christoffel symbols to obtain 

ff'a/iix, X) = -^^ig''" {^aqpvix, x) + ^ pQavix, X) - Vyq^fi(x, x)) + 0{K^) , (3.64) 

and the covariant derivatives are defined with the Riemannian metric. 
Thus, the anomalous acceleration is 

A^nom = ^^g'" (2V,qpAx, X) - Vyq^p{x, x)) xTxP + 0(^2) , (3.65) 

3.2 Uniform Model: Second Order in k 

So, in this section we will simply assume that all mass parameters are equal, that 
is, 

m 

m; = — . (3.66) 

In this case the Finsler function F{x,^) is given by (12.351) . By using the decom- 
position of the matrix- valued metric and the fact that trh^^ = we get the Finsler 
function 

Fix,e = VH?F(i - k'Is'^^^^^x)^^^ + 0{k') , (3.67) 



17 



G^"'(x, = g^^ix) - K^-S^^^^xY-^ + 0(x^) , (3.70) 



G^y(x, U) = g^y(x) + ^^J^^yaysWl^r + ^^^^^ " ^'^'^^^ 



and the Hamiltonian 

H(x,0 = ^\^\'\l-x'^S''^'"'(x)^^^^^j + 0(x'), (3.68) 
where 

S^"""^ = ^trih^^h"^). (3.69) 
By using the above, we compute the Finsler metric 

'\s^^^^x)^Y^' 

and its inverse 

1 m"m^ 
We also note that 

1^1' = + >c^:^S,y,„(x)—^ + 0(^3) . (3.72) 
Thus, our Lagrangian is 

F(x,^(x,x)) = VHiF|l + x"^5^v^;;(x) -^''-^y^ j + 0(x') , (3.73) 
Finally, we compute the ChristofFel symbols to obtain 

ff'afiix, X) = X" ^aSpypa + - V,S ^ + 0(x^) . (3.74) 

Thus, the anomalous acceleration is 

x^ / \ xFxTX^xP 

A^anom = -yg^" (2V„5^vpcr " yS aPpa) + 0(;^3) , (3.75) 

Notice that with our choice of the parameter t we have F{x,^{x,^)) = 1, and, 
therefore, in the equations of motion we can substitute with the same accuracy 

= -1 + o(k^) , \x\- = -1 + 0(x^) . (3.76) 

Therefore, we obtain finally 

2 

A'^anom = -jg""' (2V,5^,p, - V,S„^p,) x^x'^x'^x^ + 0(x') . (3.77) 
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3.3 Non-commutative Corrections to Newton's Law 

Now, we will derive the non-commutative corrections to the Newton's Law. We 
label the coordinates as 

/ = t, =r, x^ = e, x^ =(p, (3.78) 

and consider the static spherically symmetric (Schwarzschild) metric 

ds'^ = -U(r)df- + U~\r)dr^ + r^{de^ + %m^ed(i) , (3.79) 

where 

U(r) = 1 - - , r„ = 2GM , (3.80) 
r 

and M is the mass of the central body. It is worth recalling that here t is the 
coordinate time. In the previous sections we used t to denote an affine parameter 
of the trajectory that we agreed to choose to be the proper time. In the present 
section we use r to denote the proper time and t to denote the coordinate time. 

The motion of test particles in Schwarzschild geometry is very well studied in 
General Relativity, see, for example IfTOl . Assuming that the particle moves in the 
equatorial plane 6 = n/2 away from the center, that is, dr/dr > 0, the equations 
of motion have the following integrals [fTOll 

,0 = ^ = ^J_, (3.81) 
dr mU (r) 

x' = f^ = ^l^-\l + ^-\U(r), (3.82) 

de 

dr 




. de 

x^ = — = 0, e=-, (3.83) 



.,^d^ = (3.84) 
dr m 

(3.85) 

where m, L, and E are the mass of the particle, its orbital momentum and the 
energy. 

In the non-relativistic limit for weak gravitational fields, assuming 



E = m + E' , 



(3.86) 
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with E' « m, and r » one can identify the coordinate time with the proper 
time, so that 

i° = — = 1 . (3.87) 
dr 

Further, for the non-relativistic motion we have r, r6, r(p « 1, and the radial 
velocity reduces, of course, to the standard Newtonian expression 



dr I IE' L2 1 r, 

dr y m r 



x' = — = A 7^ + -' (3-88) 



which for L = becomes 



^, * /2£:^^^ (3.89) 

dr y m r 

It is worth stressing that the anomalous acceleration due to non-commutativity 
in the non-relativistic limit can be interpreted as a correction to the Newton's Law. 
Assuming that a particle is moving in the equatorial plane, 6 = n/2, with zero 
orbital momentum, (p = const, the equation of motion is 

d^r d 

'^^+A'-,„om, (3.90) 



r2 



where in the uniform model 



A\,,^ = -drS"'^'' + 0{k'), (3.91) 



with 5°"°° = j^ixh^^h^^, and in the non-uniform model 

A^anom = -^5.^°° + 0(x2), (3-92) 

with being the component of the tensor q^^ defined by (13.591 ). This gives the 
non-commutative corrections to Newton's Law: in the uniform model, 

Kff(r) = -— -^5°«°V) + 0(;.^), (3.93) 
r 8 

and, in the nonuniform model, 

GM K no 9 

Kff(r) = + -^""(r) + Oix") . (3.94) 

r 2 

Here, of course, the tensor components S and q'^'^ should be obtained by the 
solution of the non-commutative Einstein field equations (in the perturbation the- 
ory). 
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4 Noncommutative Einstein Equations 

The dynamics of the tensor field d'^ is described by the action functional of Matrix 
Gravity. As was outlined in the Introduction, there is no unique way to construct 
such an action functional: there are at least two approaches, one ^ based on 
a non-commutative generalization of Riemannian geometry and another one ^ 
based on the spectral asymptotics of a non-Laplace type partial differential op- 
erator. The exact non-commutative Einstein equations for the action functional 
proposed in [2J were found in our recent work |[6l . The equations of motion for 
the spectral approach were obtained within the perturbation theory in our paper 
II3. In the present work we are using the approach of 121. We should mention that 
in the perturbation theory the difference between these two approaches consists in 
just some numerical parameters of the action; the general structure of the terms is 
the same. 

In the following we will give a very brief overview of the general formalism, 
more details can be found in |I21[6l. We define the matrix- valued tensor b^^y by 

fl^^Z7,, = 5^I, (4.1) 

the matrix- valued connection coefficients ^^ap by 

^"a, = hbAaia^^dyaP" - aP^dyd"" - a''''dya"P)bp, , (4.2) 
and the matrix- valued Riemann curvature tensor 

^ afiv — d^^^ ay dy^^ a/i f}^^^^ av »^ fiv^^afi • (4.3) 

Next, we define a matrix-valued density 

p= f ^cxp[-a''%^y]. (4.4) 

This enables us to define the action of Matrix Gravity as follows 

^MG(«) = / Re ifrp (a^^T^Vv - 2A) , (4.5) 

M 

where G is the Newtonian gravitational constant and A is the cosmological con- 
stant. 
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Of course, in the commutative limit all these constructions become the stan- 
dard geometric background of General Relativity. The tensors and bf^y be- 
come the contravariant and covariant Riemannian metrics, the coefficients Jl^^iv 
become the Christoffel symbols, the tensor 'R^^f^y becomes the standard Riemann 
tensor and the action of Matrix Gravity becomes the Einstein action functional. In 
the presence of matter one should add to this functional the action of the matter 
fields and particles. The non-commutative Einstein equations were obtained in 
|6J. These equations, in full generality, are a complicated system of non-linear 
second-order partial diff"erential equations. Their study is just beginning. Of 
course, it would be extremely interesting to obtain some simple exact solutions. 

In the present paper we study the effects of these equations in the simplest 
possible case restricting ourselves to a commutative algebra. The commutativity 
assumption enormously simplifies the dynamical equations. In this case they look 
exactly as the Einstein equations in the vacuum 

"R,. = A^^v , (4.6) 
where !^^v is the matrix-valued Ricci tensor defined by 

f^uv = ^"uffV • (4.7) 



4.1 Static Spherically Symmetric Solutions 

We study, now, the static spherically symmetric solution of the equation (14.61 ). We 
present the matrix-valued metric a^'' by writing the "matrix-valued Hamiltonian" 



sm 9 



(4.8) 



or the "matrix- valued interval" 

b^ydxi'dx" = A'\r)dt^ + B-\r)dr^ + I (^dO^ + sin^ ed(p'^) , (4.9) 

where the coefficients A(r) and B{r) are commuting matrices that depend only on 
the radial coordinate r. This simply means that we choose the following ansatz 



^00 



-.22 



= A 



a^^ = B , 



33 

a 



1 



sin^ 6 



I . 



(4.10) 
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Next, by computing the connection coefficients ^"/^v and tfie matrix-vaiued 
I^cci tensor we obtain tlie equations of motion 



13 1 1 

-A-^A" - -A-^{A'f + -A ^A'B^B' + -A~^A' 
2 4 ' ' 4 r 



AA 



-1 



-A-' A" - -A-^{A'f + -A-^A'B-^B' - -B'^B' = Afi" 
2 4 ' ' 4 r 

--B' -B + -BA-^A' + I = Ar^ ■ I , 
2 2 

sin^ e ^22 = Ar^ sin^ 6-1. 



wtiere tfie prime denotes differentiation witli respect to r. 
By using tlie equations (14.111) and (14.121) we find 

A + B~^B' = ; 

the general solution of this equation is 

Air)B(r) = C, , 



(4.11) 
(4.12) 

(4.13) 
(4.14) 

(4.15) 
(4.16) 



where Ci is an arbitrary constant matrix from our algebra. We require that at the 
spatial infinity as r ^ 00 the matrices A and B and, therefore, the matrix C as 
well, are non-degenerate. 

By using this relation we obtain further from eqs. (|4.12l) and (|4.13l) two com- 
patible equations for the matrix B 



B" + -5' + 2A = . 



and 



rB' + B = (1 - Ar^)!. 
The general solution of the eq. (I4.18D is 



1 



5(r) = 1 - -Ar" I + -C2 , 



1 



(4.17) 
(4.18) 

(4.19) 



where C2 is another arbitrary constant matrix from our algebra. It is not difficult 
to see that this form of the matrix B also satisfies the eq. (14.171) . The matrix A is 
now obtained from the equation (14.161) 



A(r) = Ci 



1 - -Ar^ 1 1 + -C2 



1-1 



(4.20) 
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We will also require that in the limit x —> we should get the standard 
Schwarzschild solution with the cosmological constant 

B(r) = -A-\r) = |l - ^Ar^ - ^^jl, (4.21) 

where is the gravitational radius of the central body of mass M, 

rg = lGM, (4.22) 
that is, in that limit the matrices Ci and C2 should be 

Ci=-I, C2 = -rgl. (4.23) 

4.2 2x2 Matrices 

To be specific, we restrict ourselves further to real symmetric 2x2 matrices gen- 
erated by 

In this case the constant matrices Ci and C2 can be expressed in terms of four real 
parameters 

Ci=aI + eT, C2=t£ + LT, (4.25) 

where 9 = x9 and L = xL are the parameters of first order in the deformation 
parameter x. Here the parameters a and 9 are dimensionless and the parameters fx 

and L have the dimension of length. 
Then the matrix B(r) has the form 

B(r) = 1 1 - - Ar^ + ^ 1 1 + -r . (4.26) 
\ 3 r r 



Next, noting that = I, and by using the relation 

(fli + bT)-^ = (fli - br) , (4.27) 

- b'^ 

we obtain the matrix A(r) 

A(r) = ^(r)I + i/r(r)T, (4.28) 
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where 

^(r) = ^—r^ , (4.29) 

^^('') = 7^ -TT^ ■ (4.30) 

The parameters a, 6, fj. and L should be determined by the boundary conditions 
at spatial infinity. The question of boundary conditions is a subtle point since we 
do not know the physical nature of the additional degrees of freedom. We will 
simply require that the diagonal part of the metric is asymptotically De Sitter. 
This immediately gives 

a = -l. (4.31) 

Now, we introduce a new parameter 

ro = |Ar'^', (4.32) 

and require that for rg « r « r^, the diagonal part of the metric, more precisely, 
the function ip(r) is asymptotically Schwarzschild, that is, 

^(^) = _1 _ ^ + of^] + of^j . (4.33) 



This fixes the parameter // 

H = -rg + 9L. (4.34) 

The parameters 6 and L remain undetermined. 
Finally, by introducing new parameters 

p = (1 + ^)L - 9rg (4.35) 

n = rg-(e± 1)L (4.36) 
we can rewrite our solution in the form 

-r(r- \Ar^ -r„ + 20L) 

fi^) = r \ u . , 1 ' (4.37) 

r - ^ Ar^ - r_ r - 5 Ar^ - r+ 

r\e(r- {Ar^]+p] 

Hr) = 7 — T • (4-38) 

[r - iAr3 - r_J [r - |Ar3 - r+J 
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Of course, as x ^ both parameters L = xL and 6 = x9 vanish and we get the 
standard Schwarzschild solution with the cosmological constant. 

Notice that the matrix- valued metric a'^^ becomes singular when the matrices 
A and B are not invertible, that is, when 

detA(r) = 0. (4.39) 

The solutions of this equation are the roots of the cubic polynomials 

r - ^aP - r- = and r - ^Ar^ - r+ = (4.40) 

Recall that the standard Schwarzschild coordinate singularity, which deter- 
mines the position of the event horizon, is located at r = r^. The presence of 
singularities depends on the values of the parameters. We analyze, now, the first 
eq. in (14.401) . In the case A < the polynomial has one root if r_ > and does not 
have any roots if < 0. In the case A > it is easy to see that: i) if r_ > (2/3)ro, 
then there are no roots, ii) if < < (2/3)ro, then the polynomial has two roots, 
and ii) if r_ < 0, then the polynomial has one root. The same applies to the second 
eq. in (14.401) . 

We emphasize that there are two cases without any singularities at any finite 
value of r. This happens if either: a) A < and r+ < 0, or b) A > and 
r+ > (2/3)ro. This can certainly happen for large values of \6\ and \L\. In particular, 
if 6 and L have the same signs and 

1^1 > 1 + ^ , (4.41) 

then both r+ are negative, r+ < 0, and if 6 and L have opposite signs and 

2 

|e|>l + ^^^|^, (4.42) 

then r+ > (2/3)ro. This is a very interesting phenomenon which is entirely new 
and due to the additional degrees of freedom. 

We would like to clarify some points. The parameters /i, introduced in the 
previous sections describe the properties of the test particle, that is, the matter. The 
parameters 6 and p introduced in the static and spherically symmetric solution of 
non-commutative Einstein equations describe the properties of the gravitational 
field, that is, the properties of the source of the gravitational field, that is, the 
central body. The parameters 6 and p are not related to the parameters yU;. 
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5 Anomalous Acceleration 

In this section we are going to evaluate the anomalous acceleration of non-relativistic 
test particles in the static spherically symmetric gravitational field of a massive 
central body. 

All we have to do is to evaluate the components of the anomalous acceleration 
(|3.77l) . As we will see the only essential component of the anomalous acceleration 
is the radial one A' anom- AH other components of the anomalous acceleration are 
negligible in this limit. As we will see below, the anomalous acceleration is caused 
by the radial gradient of the component of the matrix-valued metric, which is 

xh^ = tl/(r)T, (5.1) 

where i/^(r) is given by (14.381 ). Our analysis is restricted to the perturbation theory 
in the deformation parameter x (first order in x in the non-uniform model and 
second order in x in the uniform model). That is, we should expand our result 
in powers of p and 6 and keep only linear terms in the non-uniform model and 
quadratic terms in the uniform model. 

For future use we write the function i/f{r) in the first order in the parameter x 

«A(r) = ^ + 0(x^) , (5.2) 

(r - iAr3 - r,) 



and for r « ro 

r(9r + p) ^ 
'A('-) = T ^ + 0(x'), (5.3) 



and, finally, for rg « r « tq. 



ilr(r) = e+- + 0{x'-) . (5.4) 
r 

We would like to emphasize at this point that the perturbation theory we are 
going to perform is only valid for small corrections. When the corrections become 
large we need to consider the exact equations of motion (|2.51l) . 

5.1 Uniform Model 

In the non-relativistic limit the formula for the anomalous radial acceleration 
(1MB gives 

A\nom = ^iA(r)(A'(r) + 0(^3) . (5.5) 
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The derivative of the function ij/{r) is easily computed 

= co{r)ilj{r) , 

where 

1 - Ar2 



1 0(1 -Ar2) 
ti;(r) = - + 



1 -Ar2 



e(r-|Ar3)+p r-iAr3-r_ r - ^Ar^ - r+ 
Thus, we obtain finally 

1 



(5.6) 



(5.7) 



A'' 



-r(r)a>(r) + C>(%0 



(5.8) 



Recall that the parameters p and 6 are of first order in x. Strictly speaking we 
should expand this formula in p and 8 keeping only quadratic terms; we get 



1 [e[r- \Ar^)+p]r 
-2r(l - Ar^) 



r Ar - 



2r--ArM+p 



elr-^Ar^l+p 



(5.9) 



For r << ro (that is, |A|r^ << 1) this becomes 

ir(er + p) [{p + 2erg)r + pr^ - f ^A/] 



4'' 



4 (r - rg)5 



+ 6>(;^^). (5.10) 



We need to keep the term linear in A since we do not know the values of the 
parameters 6 and p. Finally, for Vg « r « ro we obtain 



-eXr\ + 0{K'). 



(5.11) 



5.2 Non-uniform Model 

Similarly, in the non-uniform model the anomalous acceleration is given by eq. 
(13.921) . In the 2x2 matrix case considered above the eigenvalues of the matrix 

h^'^^l,^, are 

^i,2 = ±^tr(/i'^V)^,^,. (5.12) 
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Therefore, 

where 
Thus 

So, we obtain 
Thus 



P(x, ^) = //I Ai + i^2A2 = r-tr (h'''T)^^^y , 



(5.13) 
(5.14) 
(5.15) 
(5.16) 



(5.17) 



Now, we recall that p and 9 are of first order in x and expand in powers of p 
and keeping only linear terms 



r 



(r - |Ar3 - r J 



3 



e\2r- ^Ar'\+p 



-2r(l - Ar^) 



9^- ^^'"^^j + P 



(5.18) 



In the case r « tq (when |A|r^ « 1) this takes the form 



anom 2^" 



[(p + 29rg)r + pVg - fMr^] 



+ 0(piO • 



(5.19) 



Finally, for r„ « r « ro we obtain 



anom 2 ' 



+ 0(?^2) , 



(5.20) 
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6 Conclusions 

In this paper we described the kinematics of test particles in the framework of 
a recently developed modified theory of gravitation, called Matrix Gravity [[H 
mill. We outlined the motivation for this theory, which is a non-commutative 
deformation of General Relativity. Matrix Gravity can be interpreted in terms of a 
collection of Finsler geometries on the spacetime manifold rather than in terms of 
Riemannian geometry. This leads, in particular, to a new phenomenon of splitting 
of Riemannian geodesies into a system of trajectories (Finsler geodesies) close 
to the Riemannian geodesic. More precisely, instead of one Riemannian metric 
we have several Finsler metrics and different mass parameters which describe 
the tendency to follow a particular Finsler geodesies determined by a particular 
Finsler metric. As a result the test particles exhibit a non-geodesic motion which 
can be interpreted in terms of an anomalous acceleration. 

By using a commutative algebra we found a static spherically symmetric so- 
lution of the modified Einstein equations. In this case a completely new feature 
appears due to the presence of additional degrees of freedom. The coordinate 
singularities of our model depend of additional parameters (constants of integra- 
tion). Interestingly, there is a range of values for these free parameters in which no 
singularity occurs. This is just one of the intriguing differences between Matrix 
Gravity and General Relativity. 

The description of matter in Matrix Gravity needs additional study. In this pa- 
per we studied just the behavior of classical test particles. We propose to describe 
a gravitating particle by several mass parameters rather than one parameter as in 
General Relativity. We considered two models of matter: a uniform one, in which 
all mass parameters are equal, and a non-uniform one, in which the mass param- 
eters are different. The choice of one model over the other should be dictated by 
physical reasons. It is worth emphasizing that in the generic non-uniform model 
the equivalence principle is violated. 

The interesting question whether the matter is described by only one mass 
parameter or more than one mass parameters as well as the more general question 
of the physical origin of multiple mass parameters requires further study. Since 
we do not know much about the physical origin of these mass parameters masses, 
we do not have to assume that they are positive. We do not exclude the possibility 
that some of the mass parameters can be negative or zero. This would imply, of 
course, that in this theory there is also gravitational repulsion (antigravity). This 
could help solve the problem of the gravitational collapse in General Relativity, 
which is caused by the infinite gravitational attraction. 
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The next step of our analysis of the phenomenological consequences of Matrix 

Gravity is to apply the kinematic model developed in the previous sections to the 
study of such effects as the motion of Pioneer spacecrafts (Pioneer anomaly) and 
galactic rotations (dark matter). It would be very interesting to understand if the 
anomalous acceleration of the spacecrafts and the flat rotation curves of galaxies 
can be explained without the concept of dark matter. 
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